Macroscopic signature of protected spins in a dense frustrated magnet 
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The inability of systems of interacting objects to satisfy all constraints simultaneously leads to 
frustration. A particularly important consequence of frustration is the ability to access certain 
protected parts of a system without disturbing the others. For magnets such "protectorates" have 
been inferred from theory and from neutron scattering, but their practical consequences have been 
unclear. We show that a magnetic analogue of optical hole-burning can address these protected spin 
clusters in a well-known, geometrically-frustrated Heisenberg system, Gadolinium Gallium Garnet. 
Our measurements additionally provide a resolution of a famous discrepancy between the bulk 
magnetometry and neutron diffraction results for this magnetic compound. 

PACS numbers: 75.40.Gb, 75.45.+j, 75.50.Dd, 75.50.Lk 



One important property of frustrated systems is that 
the number of most probable (lowest energy) configura- 
tions is not finite and grows exponentially with the sys- 
tem size, resulting in finite zero-temperature ( T) entropy. 
A number of real magnets, which are among the simplest 
physical realizations of frustrated systems exhibit non- 
zero entropy [ESQ as T^O. What has been missing 
is a simply measured, macroscopic, dynamical effect as- 
sociated directly with decoupled spin subsets. We report 
here the discovery of such an effect in a dense, ordered 
network of spins. Our result shows that magnetic de- 
grees of freedom can be 'protected' [H from each other 
in a dense magnetic medium. It also provides a sharp 
distinction between a well-known magnet, Gadolinium 
Gallium Garnet, where the incompatible constraints are 
due to lattice geometry, and conventional spin glasses, 
where they are due to disorder. 

The archetypical frustrated magnet is the triangular 
antiferromagnet (AFM) (Fig. la) consisting of Ising 
spins. If, on a particular triangle, a spin is antiparal- 
lel to its two neighbors, up and down orientations have 
equal energy. This spin is a protected degree of free- 
dom because it can flip as if it were decoupled from its 
environment. Such degrees of freedom involving larger 
groupings of spins can be found for more complicated lat- 
tices (Fig. ib) HUBS. Celebrated three-dimensional 

generalizations of the Kagome lattice include pyrochlorcs 
|9| and garnets. Here we focus on Gadolinium Gallium 
Garnet (Gd3Ga50i2 or GGG), a cubic compound con- 
taining magnetic Gd ions on corner-sharing triangles that 
form two interpenetrating networks (Fig. lc). The Weiss 
temperature ©weissi proportional to a combination J of 
the exchange (1.5 K) and dipolar (~ 0.7 K) couplings 
between spins, is approximately 2 K [13] .The single ion 
anisotropy is less than 0.040 K, making this a Heisenberg 
system. The geometric frustration postpones magnetic 
order to T <0.180 K, where the bulk magnetic response 
111 ] begins to resemble that of a magnetic glass. How- 




FIG. 1: (a) Spins at vertices of a triangular Ising antiferro- 
magnetic (AFM) lattice (labeled "?") are geometrically frus- 
trated (b) Heisenberg AFM on a Kagome lattice. The "local" 
degree of freedom corresponds to a collective rotation of the 
red and blue spins around the black spin, (c) One of the 
Gd sublattices of GGG, showing the ten-membered ring of 
spins which can rotate about the black spin direction, form- 
ing a local degree of freedom analogous to that in (b). (d-f) 
Schematics of spin triads to understand response of local de- 
grees of freedom to small external fields h ac . The black spin 
sublattice is assumed to be rigid, (a) and (b) are adapted 
from 0]. 



ever, muon spin relaxation does not display the behavior 
conventionally associated with spin glasses below T g [13] • 



Additionally, neutron scattering 13| reveals sharp mag- 
netic diffraction peaks superposed on a spin liquid-like 
structure factor. This, together with recent theory 14j . 
suggests that in the pure limit, conventional long range 
AFM order, rather than a disordered spin glass, coex- 
ists with the spin fluid that arises from local degrees of 
freedom as illustrated in Fig. lc. 
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FIG. 2: (a) Phase diagram (following [TJl). Our study focuses 
on the spin liquid phase with H <4 kG< 6 kG. The red 
dashed lines indicate the temperatures along which we have 
conducted detailed dc (H) studies, (b) X "t as a function of 
T and / at H — in the linear regime, (c) H dependence of 
x" ■ The low / plateau in x"(f) connected to spin ordering at 
0.100 K disappears rapidly with the application of H. 



To probe the true nature of the ground state (Fig. 2a) 
we measured the dynamic magnetization M for oscillat- 
ing drive fields h ac cos(2irft) down to T=0.050 K as a 
function of the amplitude h ac and frequency /. The sam- 
ple is a single crystal of GGG (long axis [100]) grown 
by the Czochralski method and obtained from Princeton 
Scientific. In the limit of small h aci M is a linear function 
of h ac and the magnetic susceptibility x(/)=^MaC In 
fig. 2(b) we plot the imaginary part x"(f) against / and 
T for /i ac =0.04 G, which we have established (see below) 
to be in the linear regime. At the high temperature end 
(0.500 K, black), X "(f) has a peak at 80 Hz with a FWHM 
of ~ 1.3 decades in /, implying that the dynamic response 
is described by a distribution of relaxation times [15j . As 
T is lowered, the peak height of x"(f) increases while 



the basic shape and f pea k remain almost the same, until 
0.110 K, where the shape alters drastically. There is no 
longer a peak in \"{f), which plateaus at low /. The flat, 
/-independent tail for /— *0 originates from 1// noise in 
M, which is characteristic of the dynamical scale invari- 
ance characteristic of second order phase transitions as 
well as certain ordered phases with complex hierarchies 
of ground states, such as spin glasses [lj|. Cooling be- 
low 0.100 K reduces x"(f) and its shape returns to that 
seen above 0.110 K, again with f pe ak—8Q Hz. This is 
neither expected nor seen for conventional spin glasses 
[l6j |. What we observe looks more like a transition to a 
conventional ordered state, where there is critical slow- 
ing down and condensation of slow modes at the phase 
transition, below which the excitations harden as the or- 
der grows. The picture given by the linear susceptibil- 
ity is consistent with neutron scattering data [l3| , where 
relatively sharp magnetic diffraction peaks correspond to 
partial AFM order for T<T N = 0.100 K< 0.140 K. How- 
ever, our bulk measurements reveal something unusual, 
that the characteristic frequency below which a scale- 
invariant response is observed corresponds to an energy 
0.3 peV<fc B 7V ~ 10/^eV< J-200 ^eV. For an ideal soft 
mode transition, all of these energy scales would be of 
the same order. However, soft modes are often accom- 
panied by a 'central peak' [13], typically thought to be 
a contribution of very slowly relaxing clusters of ordered 
material in a disordered matrix above Tat. 



Given the previous work [ll| on the glass-like behavior 
of GGG, it is natural to ask whether our data are in the 
linear regime. Fig. 3a-d shows measurements of M(h ac ) 
for /=29.9 Hz at a variety of temperatures. At 0.4 K we 
are clearly in the linear regime while at Ta?=0.100 K, 
there are obvious signs of a crossover from a large, low- 
h ac response to a substantially smaller high-/i ac magnetic 
response. At the same time, the imaginary part of the 
magnetization has a strong peak, which means that dissi- 
pation is reduced on going into the high-/i ac regime. The 
highly non-linear response of GGG revealed at low h ac 
in the current experiment clearly needs to be taken into 
account when analyzing magnetic susceptibility data. In- 
deed, the AFM phase transition seen in the neutron ex- 
periments and in our linear susceptibility data can be 
masked by the non-linearities which clearly matter for ac 
magnetometry with h ac =l G [ill ]. 

Similar behavior [l8} has been seen in the dilute Ising 
salt LiHo 2; Yi_ 2 ;F4, but with the crucial difference that 
there, at high h ac , M" ^0 and saturation and phase 
locking are complete. By contrast, in concentrated GGG 
there is a persistent linear (in h ac ) background. Be- 
cause we were unable to derive an analytic expression for 
M(h ac ), we have simply evaluated the low and high-/i ac 
derivatives, dM'/d/i QC and dM" /dh ac , as a function of T, 
plotted in Figs. 3e and 3f. The high h ac response shows 
a strong peak above the ordering temperature, moving 
towards the extrapolated (from high T) Curie- Weiss sus- 
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FIG. 3: (a-d) Real (M', black, left axis) and imaginary 
(M",md, right axis) parts of the magnetization as a func- 
tion of hac at different T for /=29.9 Hz. (e) The T depen- 
dence of dM'/dhac at low (/i ac =0.04 G) and high (h ac =lA 
G) drive fields. The line (blue, with solid circles) corresponds 
to the Curie-Weiss law representing the extrapolated high- 
temperature susceptibility dominated by the frustrated units 
of GGG. (f) The analogous plot for M" . (g) The spectral 
width of the hole in the susceptibility as a function of tem- 
perature. All lines are guides to the eye. 

ceptibility [10] of GGG, while the low h ac response has its 
sharp maximum at XV, with an amplitude correspond- 
ing to 10% of the Gd ions rotating freely. 

When faced with a new dynamical phenomenon in a 
system with many degrees of freedom, it is important to 
determine whether it is due to the relaxation of many 
coupled spins, or to the motion of smaller groups of de- 
coupled spins. The classic experimental method for dis- 
crimination between the two possibilities is hole burning 
fl9j |. We apply simultaneously a sinusoidal pump, with 
amplitude 0.8 G, sufficient to approach the crossover to 
the high h ac regime in Fig. 3a-d, and a small ampli- 
tude probe (0.04 G). Fig. 4a shows how at 0.150 K, 
the pump carves out a hole in centered at f pU mp 

and with FWHM -3 Hz. No other part of %"(/) is simi- 
larly affected although at / < 250 Hz the entire spectrum 
is attenuated relative to that observed without pump- 
ing. The holes contain about 0.5% of the total spectral 
weight seen at 0.150 K, a macroscopic number ~ 10 17 
spins cm -3 . We can burn similar holes at any frequency 
/ < f peak . At 0.095 K (Fig. 4(b)), just below T N , we see 
a sharper, deeper hole, which then becomes broader and 
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FIG. 4: (a-c) \" as a function of / at different T bracketing 
the ordering temperature in the linear (black) and non-linear 
(blue) regimes. The non-linear susceptibility is measured with 
a pump-probe combination. /j,« m j,=29.9 Hz (blue arrow), (d) 
The holes in (a-c) magnified, sharpest at ordering, (e) Hole- 
burning at 0.100 K at nonzero H. The ability to burn holes 
disappears by _ff=4 kG (data shown in the bottom panel of 
(e) with an added offset for clarity). 



shallower well within the ordered phase, at 0.050 K. Fig. 
3(g) shows how the spectral hole width, inversely propor- 
tional to the quality Q of the oscillators responsible for 
the excitations at 29.9 Hz, has a sharp minimum at Tjy 
We have also applied a dc bias field H to approach the 
different AFM phase (see Fig. 2c) revealed previously 
above H—6 kG. As shown in Fig. 4e, the sharp hole is 
robust up to at least 0.5 kG, but it is greatly broadened 
and suppressed at 2 kG and unrecognizable at 4 kG. 

What could give rise to the apparent high-Q oscilla- 
tors in a dense frustrated magnet such as GGG? Most 
likely are the red and blue spin clusters, shown in Fig. 
lc, which can be rotated at no cost relative to an ordered, 
surrounding backbone of black spins. A small external 
field, such as h ac , will break this rotational symmetry 
as long as it is not parallel to the field produced by the 
backbone. In particular, once the anisotropy field Ha is 
exceeded, there will be a spin flop, shown schematically 
for the spin triad in Fig. Id, with the red and blue spins 
at an angle of 30° relative to the black spin and simulta- 
neously orthogonal to h ac - There will then be some small 
additional rotation (89) of the red and blue spins towards 
each other and h ac . In the pure limit, where there are no 
random internal fields, there will a much smaller response 
(Fig. le) for h ac < Ha than for h ac > Ha- Indeed, the 
imaginary part of the bulk response should be rigorously 
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zero [20[ for h ac — > for a clean Heisenberg magnet (with 
Ha = 0) consisting of AFM triangles. 

Disorder changes the situation markedly and the re- 
sponse can be stronger for h ac < Ha than for h ac > Ha- 
For example, if there is a missing magnetic atom on the 
ten-membered ring of Fig. lc, the ring will have a net, 
uncanceled moment M u which is orthogonal to the black 
spins, and before the spin flop transition occurs the main 
effect of small h ac will be to orient M u , a situation 
which can be understood by simply removing one spin 
from Fig. le to obtain Fig. If. Above h,A, the spin 
flop transition will still occur and M u will be perpendic- 
ular to h ac , but with a correspondingly reduced magnetic 
response dM/dh ac . This leads to precisely what we ob- 
serve in Fig. 3. Dissipation is reduced at the crossover 
between the spin-flop regimes because continuous rota- 
tion of moments orthogonal to h ac , as illustrated in Fig. 
Id, is intrinsically barrier-free, in contrast to the flipping 
of moments for h ac < Ka as the applied field oscillates. 

We can estimate the anisotropy field Ha to be close 
to the value of h ac where the susceptibility crosses be- 
tween the linear and non-linear regimes and from Fig. 3 
this is approximately 0.8 G. Yia also leads to an under- 
standing of the observed characteristic frequency f pea k 
which corresponds to the most probable relaxation rate, 
typically of order J, translating to 50 GHz for GGG. 
Because we are dealing with the relaxation of groups 
of spins rather than individual moments, the observed 
fpeak —80 Hz is much smaller. In this case, the relax- 
ation rate can be estimated using a WKB approxima- 
tion [2l| as f P eak=foe~ s / h , with the tunneling action 
S/h=Ns^2k/ AE for magnetic clusters with N spins 
of magnitude s, an anisotropy energy k, and an energy 
barrier AE, due e.g. to interactions with neighboring 
clusters [2l|. If we set the attempt frequency f —J and 
JV=10, the ring size of Fig. lc, we can obtain f pe ak—80 
Hz, in agreement with the experiment provided that 
2k/ AE - 1/3. If we take k =0.4 mK =gs[i B h A , then 
AE ~3 mK, which is of the order of the next nearest 
neighbor exchange interaction [Io| . 



Inevitable longer-range interactions [14j , beyond indi- 
vidual triangles, fill out the picture. In particular, at the 
lowest temperatures the fluctuating, uncompensated mo- 
ments coexist with the unsaturated AFM order revealed 
by neutrons, and via the further neighbor interactions, 
display dynamics conditioned by the static mean field 
from the background antiferromagnetism. If we impose 
static ferromagnetic correlations via an external field, we 
see the same reduction in Q that is produced by low- 
ering T below Tm for H—Q. At TV, the static mean 
field from the other clusters is by definition zero, mean- 
ing that the defect-centered clusters are truly protected 
and independent of the AFM backbone, so that here we 
obtain the highest possible Q. On further warming, the 
fluctuation rates of the backbone move through the hole 
burning frequency, coupling all modes together, and so 



Q deteriorates again. 

We have produced an experimental resolution of the fa- 
mous discrepancy between bulk magnetometry and 
magnetic neutron diffraction [13( for GGG. The transi- 
tion temperature identified previously as the entry point 
to a spin glass is actually above the condensation tem- 
perature Tjv for a continuum of ultrasoft modes in the 
low frequency, linear magnetic response. Completely con- 
trary to conventional phase transitions where modes are 
maximally coupled at a phase transition, at XV, these 
oscillators achieve maximum decoupling as measured by 
the quality factors established by magnetic hole burning. 
Our discovery that the characteristic frequency f pea k, for 
the rollover to a x" (/) decreasing with / , is independent 
of T and so is not thermally activated, indicates that the 
protected, defect-nucleated degrees of freedom are actu- 
ally behaving as quantum rather than classical objects. 
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